LPT 12-79 
LAL 12-274 



Impact of 5 — > KqI^I on the New Physics search 
in 5 ^ KH-^t decay 

^ Damir Becirevic" and Andrey TAYDUGANOV"''' 

° Laboratoire de Physique Theorique (Bat. 210) ^ 
Universite Paris Sud, Centre d'Orsay, F-91405 Orsay Cedex, France 

^ ^ Laboratoire de rAccelerateur Lineaire, Centre d'Orsay, Universite de Paris-Sud XI, 

^ B.P. 34. Batiment 200, 91898 Orsay Cedex, France 

(N 

Abstract 

We discuss the uncertainty related to the amount of unwanted B — > K^{Ktt)£~^£~ events 
^ ^ in the sample of B — >■ K* {Ktt)£~^£~ ones. Those events can increase the measured differential 

O decay rate by up to 10% in the low region, and can be a source of non-negligible uncertainty 

i-C in the full angular distribution of the B — > K* {Ktt)£'^£^ decay. Although the transverse 

asymmetries should be unaffected by the presence of the S'-wave Kt: pairs, coming from the 
^NJ scalar Kq meson, we show that in practice, their normalization might be sensitive to those 

J>" events and could entail a sizable uncertainty in transverse asymmetries around = 2 GeV^. 

For other q^'s that error is under about 10%. 
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1 Introduction 

The angular spectra of B ^ K*{-^ Ktt)£'^(.~ decay contain a number of interesting observables 
which can be measured and compared with theoretical predictions made in the Standard Model 
(SM), and hopefully result in valuable hints about physics beyond the Standard Model (BSM). ^ 
To achieve that goal the main obstacle is our inability to reliably compute the hadronic quantities 
to desired accuracy. In such a situation one selects the observables that are the least sensitive 
to hadronic uncertainties that are at the same time potentially sensitive to the new physics 
(NP) signals. It appears that three transverse asymmetries, that one can build up from the q^- 
dependent coefficient functions in the analysis of the angular distribution of — )■ K* (— t- Kii)£^£~' , 
satisfy these requirements at low g^'s and are supposed to be well measured at LHCb as well as 
at Super-B and Super-KEKB/Belle II (see also ref. [2]). 



^Laboratoire de Physique Theorique est une unite mixte de recherche du CNRS, UMR 8627. 

^Recent research on the potential of angular distribution of the B — ^ K*£'^£~ decay in the search for signals of 
physics BSM can be found in ref. [1] and references therein. The work in this direction was initiated by the authors 
of ref. [3]. 
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The asymmetry A^\q'^) has been first introduced in ref. [4] and by now it is in the catalog 
of standard quantities that are expected to probe the effects of physics BSM at low energies. Its 
most significant feature is that 

so that its non-zero value would suggest that Cj ^ 0, which could only be attributed to physics 
BSM. Away from = this quantity is sensitive to other sources of NP and studying its shape at 
low g^'s could be very beneficial for either seeing the effects BSM or constraining the NP models. 

The asymmetry A'ij^\q'^) has been introduced in ref. [5, 6] and it is particularly sensitive to 
the NP phases. In particular its intercept 

hmA^ W )- |C^|2 + |c/|2' (2j 

and away from = it is sensitive to the phases coming from operators other then the electro- 
magnetic penguin operator. 

The asymmetry A^^\q'^) has been introduced in ref. [6] and is proportional to the usual 
forward-backward asymmetry Afb{q'^), except that it is not divided by the differential decay 
rate, but only by its transverse part for which the hadronic uncertainties are better controlled. It 
has the same zeros as Afb{q'^), and a well established shape in the SM. In particular the position 
of its low q^ extremum is well defined in terms of Wilson coefficients. More specifically, 



for ll = -'~, A^^\ql)=0, 

for dt. = ^TT^, (dAPiq^)/dq^ 



(3) 



ycxtr 



SO that from the position of its non-trivial zero and extremum at low q^ , one can already test the 
SM prediction ^ 

= (4) 

The shape of A^^\q^) is particularly sensitive to c'^\ 

Most of the observed events so far are B — t- K*fi^fj,~. As we are interested in exploring 
the low q'^ region it is important to include the lepton mass effects, as to properly combine 
B — )• K*fi'^iJ,~ with B — )• K*e~^e~ events. A problem that is often ignored in the literature is the 
contamination of the angular distribution of i3 — t- K*{-^ K'k)I!.'^1^ by the events coming from 
B —7- Kq{-^ KTT)i^i~, where Kq stands for a broad scalar meson resonance. This effect was 
recently studied in the experimental analysis of e.g. D — )• K*fiu decay [7, 8, 9], and was shown to 
be important. For the above mentioned asymmetries this is not a problem because the product of 
the Kq — )• Ktt decay is in its S'-wave and cannot make any impact on the B — )• K*{-^ K7r)i~^£~ 
transverse amplitudes. However, in the extraction of transverse asymmetries from the full angular 
distribution the unwanted {Ktt)s originating from B — )• KQi~^i~ are troublesome and result in 
an error that is g^-dependent and can be uncomfortably large, as we show in the following. 

In Sec. 2 we derive expressions for the full distribution of the combined B — )• K* (— )• K'ir)i~^£~ 
and B — )■ Kq{-^ K'ir)£^£~ decays, and in Sec. 3 we discuss the phenomenological consequences. 



is a ratio of form factors that will be defined below, which at low q'^'s behave like a constant. More specifically, 
(Ai(g2)/T2(g2))/(ms - mK*) ^ {V{q'^)/Ti{q^))/ {niB + niK-) G (0.17,0.23) GeV-^ See discussion in ref. [6]. 
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2 B ^ K*{^ K7t)£+£- and B K^{-^ K7t)£+£- 

2.1 Operator basis and the hadronic matrix elements 

The effective Hamiltonian used to describe the b — )• s£~^£~ decay reads [10], 

6 



-Hes 



AGf 



i=l 1=7,8,9,10 ^ ^ 



(5) 



where the twice Cabibbo suppressed contributions (oc VubV*^) have been neglected and the oper- 
ator basis in which the Wilson coefficients have been computed is [11, 12]: 



07 = 


e 
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09 = 












Oio = 









O' 



;{si^,PBb){hV), 



(6) 



with Pl,r = (1 =F 75)/2- The explicit expressions for Oi-e and Og^ can be found in ref. [11], 
together with the Wilson coefficients Ci{iJ,) which carry the SM information on physics at short 
distances. The Wilson coefficients multiplying the same hadronic matrix element are combined 
into effective coefficients [13] and in what follows whenever we write Cj^g^io, the effective coeffi- 
cients C'^ g -j^Q should be understood. Note also that in the SM the Wilson coefficients Cy_iq = 0. 

The hadronic matrix elements, obtained by sandwiching the above operators between B and 
K* and/or between B and Kq, describe the long distance physics and are parameterized in terms 
of the form factors. There are seven form factors parameterizing the B — t- K* transition matrix 
elements, namely. 



{K*ik,e)\sj^il-j5)b\B{p)) 



mB + rriK* 



+ iq>.{e* ■ q) ^ [A^{q^) - Ao{q^ , (7) 



"T-B + rriK 

where the partial conservation of the axial current provides the following relation. 



Az{q') 



2- rUB + niK* . , 2-^ ITT-B - TUk* . , 2\ 



2mK* 



Ai{q' 



2mK* 



A2{q' 



(8) 



that cancels the divergence at q^ = 0. Other three form factors parameterize the matrix element 
of the electromagnetic penguin operator, 



{K*{k,e)\sa^,q''{l + j5)b\B{p)) 



2ie^,p,e*''pPk''Ti{q^ 

,2 



+ [^li'm-B 



{e*-q) ip+k)jT2{q' 



9m 



m 



B 



m 



{p + k)^ 



K* 



Uq' 



(9) 



with Ti(0) = T2(0), ensuring that only one form factor describes the physical B — )• K*j decay. 
Concerning the decay to the scalar meson, the analogous matrix elements involve three new form 
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factors, namely, 

{K*{k)\s^,^,b\B{p)) 



-9m 



fAq' 



rriB + rriK* 



(10) 



2.2 B K*i 



and B ^ K, 



amplitudes 



With the above definitions in hands we can now write the B ^ K 



decay amplitudes as: 



M 



L,R 



{c'-^' T C- 
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mB + rriK* 



+ C. 



(+) 2m;, 



M 



L,R 



ms - rriK* 



a 



T2{q' 



M 



L,R _ Niy/N]^ I (_) (_) 

lmK*y q 



(-) 



imj, + 3mj,.-q')T2{q 



i^ml-ml.-q^){mB+mKMi{q^) 



mB + rriK* 



2 2 



-T,{q' 



M 



L,R NWXk*Nk' 



{ct''TC[,^)Ao{q'), 



where, for shortness, we write c!f^\Q = Cj^g^io ± Cy g ;^q, and ^ 



.Gp Oemy^ -r^* J. J 



2567r3 



m 



B 



A(m^,mi ,^^) = [g^ - {mB + mK,jV - {^b - mx^J'^] 



(11) 



(12) 



The above amplitudes agree with ref. [14]. As for those describing B — )• KqI'^I , they are given 
by: 



M'.'^'iq'-^ 



{ctWcto^)Uiq') + ci-hm,^^ 

mB + rriK' 



{ci-'^^C^^ml-ml-AUq^), 



(13) 



where an extra "prime" is used to distinguish them from the B — )• K* amplitudes which also agree 
with ref. [14]. The superscript L,R in the above expressions refers to the lepton pair chirality. 



^The above amplitudes are related to j4^'|j^g t(9^)i defined in ref. [12], as: A^^'j^Q ti'f') = — i\/3/8 ^x'||^o t(^^)- 
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Figure 1: Breit-Wigner function for the Kq: the thick curve corresponds to eq. (16) with Qk ~ 0, while the 
shaded area is obtained by accounting for the K-state by varying G [0,0.2]. The vertical stripe corresponds to 
ruK* ± 100 MeV. 



Finally, we should take into account the width of resonances, which we do by multiplying each 
amplitude by its corresponding Breit-Wigner function, namely 



M 



L,R I 
-L,||,0,tV 



M 



L,R I 
-L,||,0,tV 



L,R I, 2\ 



where the functions BWk*^^ ("^Ictt) given by 



BWK*{ml^)=M 



9k 



m 



K 



, - (m« - iF J2)2 - {rriK* - iTk^, /2f 



and the factor M is obtained from the normalization to unity, 



1. 



(14) 

(15) 
(16) 

(17) 



The second term in eq. (16) is the contribution of the well measured = ^0(1430), the mass 
and width of which are uik* = 1425(50) MeV and F^^* = 270(80) MeV, respectively. The first 
term in eq. (16) accounts for the -frQ(800) = k state, identified as a pole on the second Riemann 
sheet in the Kn — t- Ktt scattering amplitude, with the mass and width being = 658(13) MeV 
and Fk = 557(24) MeV, respectively [15]. To mimic its presence in the tail of the Breit-Wigner 
function for Kq, we vary the coupling < 0.2, where the upper value should be conservative 
enough as it results in a prominent second bump in BWk* ("^/Ctt) (see fig. 1). If that bump was as 
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large it would have been directly observed in experiments long before ref. [15]. We checked that 
in the region of mK-K £ [rriK* — 5,mK* + S\, with 5 ^ 100 MeV, our -BWx*,ii'*(?Ti|^^) reproduce 
the shapes of the corresponding Ktt form factors. The function (16) also reproduces well the 
findings of ref. [15], as well as those of refs. [9, 16], in the region of interest for our integration, 
rriKn G [rriK* - 5, mx* + 5] . 

2.3 Full distribution of B ^ K*{Kn)i+i- including B K*{Kn)i+i- 

Besides the usual variables, g^, 9^, 9k and cj), we now also need to consider the distribution of the 
mass of the Kn system. We have 



aq'^dmj^^a cos 9ia cos 9Kd(p 

+ [Jliq'^,mj^^, 9k) + 2J|(g^ ^i^)] cos 29 ^ + 2J^{q^,m\^, 9k) siv? 9e cos 2(P 

+2\/2J4,{q^ , mj^^, 9k) sin(26'£) coscj) + 2\f2J^{(^ ^ "i^^, 9k) sin^^ cos^ 

+2J6(g2, m|f^, 9k) cos9i + 2^/2Jj{q^ , m|f^, 9k) sin 9i sin (j) 

+2^/2 J^{q^ , m\^, 9k) sin 26*^ sin0 + 2Jg{q'^ ,mj^^, 9k) sin^ 9i sin2(/), 
where the explicit forms of the functions JiZ'^ {q^ ,'m?Kn^^K) look as follows: 

Jl^iq^ml^, 9k) =^Il2{q')\BWK^ {ml^)\'' W 9k , 



(18) 



^3,6,9('?', Ok) =—X^,Q^^{q^)\BWK' (mL)|2 sin^ 9k 



J'i{q^,m\^,9K) 



, (g^)| (mj,J|^ + 3XiV)| (m^JI^ cos^ 

+2^/3Re Zl"{q^)BWK*^{ml^)BW]^,{ni\^) cos9k\ , 



J!i{q\mi,^,9K)=^\Z'2\q')\BWK:,{mi,^)V + 2,X'2{r)\BWK^ 



+2^/3Re \z'^"{q^)BWK:, {ml^)BW],, (mD 



cos 9k \ , 



Mq'^,mK^,9K) 



2 47r\ 

Vs. 



Re 



i:iq')BWKsim'K^)BW],,imj,^) 



sm &K 



+ ^l4iq^)\BWK*{ml^)\^ sm29K 
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--(■ 

2 47r\ 

Vs. 



Re 



+ {q^)\BWK* (mL ) I ' sin 201, 



sm 



J7{q^,m\^,eK) 



2 47r \ 



Im 



+ ^X7(c?')| W;^* (mL)|2 sin 20^^ 



sin 9 k 



2 47r\ 



Im 



sin 6 k 



+ ^Z^{q^)\BWK*{ml^)\''siu2eK 



(19) 



The functions ll'^''^\q'^), expressed in terms of the amphtudes ■M.^^Q^{q'^) from eq. (11), are hsted 
in the appendix of the present paper, together with the functions x\^^' {q^) that involve the scalar 
meson contributions only (13), and the functions X^^'^"{q^) that correspond to the interference 
terms. ^ 

The transversity amplitudes obviously remain unchanged when the K-n from Kq are included 
in the decay's angular distribution. Only the scalar and the "t" amplitudes get changed. Since we 
focus on the transverse asymmetries only, we could have ignored the whole issue. However, since 
the transverse asymmetries are (implicitly) normalized to the differential B — )• K*l^i~ decay 
rate in a given g^-bin, it is important to estimate the size of the contribution of the i^'vr-pairs 
originating from Kq, relative to the dominant K* — )• Kir, inside the window around the K* mass. 

2.4 Separate distributions in 0, O^, Ok 

We now integrate the full distribution (18) to get the dependence of d?T /dq^dm\^ on the angles 
(/>, Ok separately. Above and in the following F stands for the decay width comprising both 
B^K*{^ Ktt)£+£- and B ^ K^{^ Ktt)1+1- modes. We get: 



[Sinq')-n'{e)\ \BWK',{ml,, 



dq^ dm\^ d cos 9k 3 

+ 2V3Re[(3Xf' (g2) -X|"(g2)) Wx5(m|jSW|.(mL)] cos^^ 
+3[3Xf(g2) -X2^(g2) _ 3l!{q^)+II{q')]\BWKAm%^)\^cos^eK 



(20) 



The functions lt'''H<l^) are related to the famiUar if '(g") of ref. [12] as: ^^'"^(g") = (3/8) /^'"'(g^) 
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[GeV^] 

Figure 2: Ratio between dVo/dq^, the B ^ K*{^ Kn)fi'^ii^ differential decay rate, and dV /dq^ that includes 
also the B — >■ iTg (— > KTv)fi'^ fj.~ events integrated over the interval rnKw G [niK' — 100 MeV, mjf* + 100 MeV]. We 
limit the discussion to low < rn^jj^. 



+ [Xl{q^) -X|(<?2) +2Xf(g2) - 2Zi{q')]\BWK<ml^)\'' 
+ 2XQ{q^)\BWK*{m\^)\^cose, 

+ 2[Z^,\q^)\BWK*{ml^)? + {r2{e) + 2n{q'))\BWK4^^^ (21) 



2 



+ [3X?(g2) +6Xf(g2) _ 2X|(g2)]|Wi,*(mL) 

+ ^-^Re[Xi;{q')BWKr,{ml^)BWU<.)\ cos<A 

+4 [X3('?') cos 2(/. + X9(g') sin 2</.] (mL)l' | ■ (22) 

To get the differential decay width we integrate over ah three angles and obtain 
d'^T 2 



dq'^dmj^^ 3 



|2 



+ [3Xf(<?2) -Il{q^) +6Itiq^) -2II{q')]\BWK'{ml^) 

(23) 



dq^drn]^^ dq^dm]^^ 

where we separated the part coming from the scalar resonance iVs) from the usual expression for 
B — )• K*i^£~ decay (Fq). We stress again that when integrating over mj^^ £ [(w-x* —^)'^, {riiK* + 
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5)^], besides the desired Ktt corresponding to K* , one also gets the Ktt pairs coming from Kq. 
In particular, in fig. 2 we show the ratio of dVQ/dq^ and dT/dq'^, after integrating over m'j^^ using 
5 = 100 MeV, as in experiments. We see that the inclusion of Ktt coming from Kq amounts to 
at most 10% excess with respect to the desired dTa/dq^. One should keep in mind that for this 
estimate: (i) we used the form factors from ref. [17] for which it is hard to assess uncertainties; 

(ii) we also had to use the form factors ^2,0(9^)) T^il^)-, that are prone to large uncertainties [18]; 

(iii) for the B — )■ Kq transition form factors we used the results of ref. [19], obtained from the 
standard QCD sum rule analysis of the three-point functions; (iv) we varied the parameter in 
eq. (16), that enters also eq. (23), as (7^ S [0,0.2]. Note also that the normalization condition 
in eq. (17) is an assumption that could be checked by a careful comparison of the D — t- K* 
and D — )• Kq semileptonic form factors computed in the narrow resonance approximation on the 
lattice with the experimental results presented in ref. [9]. 

3 Phenomenology 

3.1 Error on transverse asymmetries 

Expressed in terms of the functions X^^''^\q^), the three transverse asymmetries mentioned in 
introduction, are: 

4(2)f.2^ _ 4X3(g^) ,(im),„2^ _ 



4^)(g2)=— (24) 
^ ^ 3Xf(g2)-X|(g2) ^ ^ 



Throughout this paper we keep the lepton mass different from zero because the lowest bins (as 
close to as possible) are the least ambiguous as far as the identification of the potential 
NP signal is concerned, and in those bins the effect of the lepton mass is significant. Moreover, 
in order to consistently combine the samples of S — >• K*e^e~ and B — )• K* fi^ fi~ decays at low 
Q'^'s the lepton mass effect should be taken into account . 

In the massless lepton case the denominators in eq. (24) reduce to 3Xf(g^) — X|(q^) = 8X|(g^) 
and one retrieves the usual expressions quoted in e.g. refs. [12]. A common denominator to all 
three asymmetries is chosen for convenience; it consists of the transverse amplitudes only (those 
that are unaffected by the presence of Ktt pairs coming from Kq). In principle, and if the sample 
of — )• K*i^£~ is so large that all the coefficient functions I^'^^q"^) can be extracted from the 
full angular distribution (18), one can get the denominator unaffected by the presence of the 
Ktt pairs coming from Kq decay. If, instead, one proceeds by considering the distribution of the 
sample of — t- K*£~^£~ events in (p, 6i, and 6k separately, then the denominator 3X^(g^) — X|(g^) 
cannot be extracted without picking up the events coming from B — )• Kq{-^ KTT)i^i~ decay. To 
be more specific, we rewrite eq. (20) as 

(q^) + iff) cos Ok + ce^ {q^) cos^ Ok , (25) 



dq^d cos 6k 

where the coefficient functions are easily identified from eq. (20). For our purpose, the important 



^The plots presented in this work are obtained by treating the uncertainties to all the form factors and to 
parameter as uniform distributions, while those related to masses and widths of the states are treated as Gaussian. 
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Figure 3: A(g^) is a contribution to aej^{q^) in the distribution (25) of B ^ K*{^ Ktt)£'^£~ coming from the 
B -> KSi-^ Ktv)£+£- decay. 

one is 

+ [3Xf(g2)-X|(g2)] I dml^\BWK-.{ml^)? 

= [3It{q')-II{q')][l + A{q')] J dml^\BWK'{ml^)\' , (26) 

where, as before, the integration over comprises {rriK* —5)^ < '^'k-k — (''^K* +5)^. Without 
the B — )• Kq{-^ KTT)£~^i~ events the functions If^iq'^) = 0, and therefore aQj^{q^) would be the 
way to fix the denominators in eq. (24). Since this is not the case, and the functions X^2(g^) 7^ 0, 
the function A((jf^) 7^ 0. By using the same form factors as discussed above we obtain A(g^) 
shown in fig. 3, where the band of values covers the uncertainties coming from all the form 
factors, as well as from the variation of the phenomenological constant g^,, as discussed after 
eq. (17). Note that in this estimate the errors on i? — )■ K* form factors are small as only the 
transverse amplitudes are considered in this case, which do not involve ^2,0(9'^) not T-i{q^) form 
factors [6]. From the plot in fig. 3 we see that around q^ 2 GeV^ as many as 25 % events, 
recognized as ag^^^q'^) of B ^ K*£^£~ , might be coming from B — )• /Cq(— )• Ktt)£^£~ decay. There 
is nothing extraordinary about ~ 2 GeV^; it is simply a point at which the denominator of 
A(g^) has a minimum, while its numerator depends only mildly on g^. We also checked that 
the contamination of the Bg — )• (/)£~^£~ decay by the Bg — )• fo{980)£~^£~ events affects ao^{q^) 
by an excess that is similar in shape to the one shown in fig. (3) except that its maximum at 
about ~ 2.5 GeV^ is no larger than 0.15. For that estimate we used 5 = 10 MeV (as used by 
LHCb [20]), and the Bg — t- /q form factors computed in ref. [21]. 

As for the numerators in a!:^ {q^) and A^\q^) they are not plagued by — )• A'g(— )• Ktt)1'^£~ , 
and they are extracted from eq. (22) written as, 

= H^l^) + cos (t> + bl{q^) sin (t> + c^iq^) cos 2(/> + c'^{q^) sin 2ct) , (27) 
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Figure 4: Dependence of the interference terms on the Ktv mass around rriK* for — 1 GeV^ (blue), and 
= 5 GeV^ (red). Left plot corresponds to bxiq^ ,m^^) which, when integrated over m|f^, leads to bK{<l^) in 
eq. (25). Similarly, in the right plot we show 6^(g^,m|f^) that, after the integration in m|f^ between (mx* — S)^ 
and (tUK' + (5)^, gives the function b^(g^) used in eq. (27). For convenience, both functions are normalized to the 
total decay width. 



after identifying 4X3(^2) = 37rc^(g2), and 4X9(^2) = 37rc^(g2). Note also that 2Tra^{q'^) = dV/dq^. 
Concerning the numerator of A^^\q^), it is extracted in the usual way from eq. (21), 

-To-, T = (^Ot {q^) + ht (q^) cos 9e + ce, (q^) cos^ , (28) 

dq^dcosOi 

withX6(9^) = be,/2. 

In concluding this subsection, we repeat that the B — )• Kq{-^ K7r)£'^i~ events are unavoidable 
when studying the transverse asymmetries from the angular distribution of i? — )■ K* (— )■ K7r)i~^i~ 
decay, due to the practical impossibility of extracting the denominators in eq. (24) without 
getting any contributions coming from the functions If'{q'^). Any other choice of normalizing 
the transverse amplitudes in eq. (24) would be as much (if not more) affected by the B — )• Kq (— )• 
KTT)i~^i~ events. The only way to circumvent the problem discussed above is to consider the full 
distribution (18) and extract the coefficient functions I^''^{q'^) from the simultaneous fit of the 
angular dependence in 9k, and (f). In doing so one would have to work in the massless lepton 
limit, which means for q^ > \ GeV^. 



3.2 The interference terms 

It is interesting to comment on the terms arising from interference between the B — )• K*{-^ 
K'K)i~^£~ and the B — )■ Kq{^ Kir)£~^£~ amplitudes, that we denoted as bgj^{q'^) in eq. (25) [ex- 
plicitly given in eq. (20)], and by b^^^{q'^) in eq. (27) [and explicitly in eq. (22)]. Before integrating 

in mj^^ the functions ^SkC'?^' "^IfTr) ^^"^ ^^(9^' "^ifTr) ^ ^^^y good approximation- odd 

with respect to mj^t (c.f. fig. 4). This is a consequence of the fact that they are proportional 
to the real part of the Breit-Wigner function BWK*{m'j^^) which changes the sign at m\,. The 
contribution of BWK*{'m'j^.j^) is too small to make any significant impact when integrating over 
mj^^ e [{niK* — 5)2, (m/^. + 6)'^] with 6 ^ 100 MeV, and therefore bgj^{q'^) and ^^(g^) are not 
exactly zero but very close to it. The non- vanishing interference term is b'^(q'^), i.e. the one that 
involves the imaginary part of the Breit-Wigner function and is therefore an even function with 
respect to mj^t , as shown in fig. 4. 
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Figure 5: Dependence of the interference term b^{q'^) (divided by the positive dehnite function a^{q'^) ) on (f . It 
crosses zero at qf^^ = 1.3(3) GeV^, that coincides with the low q^ extremum of A!f'^\q^) , indicated in eq. (3). 

An interesting feature of the function b^{q'^) is that it has zero at the same point at which 
the asymmetry A^'^\q'^) reaches its extremum. To see that it suffices to write 



h%{q') 



dml^ Im [Zil{q^)BWK*^{ml^)BW*K*{ml^)\ 



oc Re 
oc Re 



Re 



{-) 



C 



10 



(_)2mb' 



oc Re M^* 



which is zero for 



fl2 

Hint 



Re a 



(-) 
7 



a 



(29) 



(30) 



as illustrated in fig. 5. Therefore a study of the interference term V^{q'^) can be informative. In 
practice, the b^{q'^) signal can be made more significant if a larger value of 6 is chosen. 

4 Summary 



In this paper we discussed the impact of — t- Kq(-^ KTT)i^i~ events on the angular distribution 
of the B — )• K* (— )• KTr)£~^i~ decay. A study of the latter is expected to lead to important clues 
about the potential signals of physics BSM. In particular the transverse asymmetries are supposed 
to be unaffected by the 5-wave Kir pairs. We show that in practice, however, the extraction 



of these asymmetries at low q s is plagued by the B — )• Kq\ 
corresponding error is under 10% for < 1 QeV^ and for 4 GeV^ < q"^ < m 



KTr)£~^i events and that the 
j^^, while it can be as 

large as 25% around ~ 2 GeV^. Similar situation occurs in the Bs — s- ipi^ K^)i^i~ decay 
except that the effect of the Bg — )■ /o(— )• K~^K~)i'^£~ decay is smaller. It remains under 15% 
around q'^ ~ 2.5 GeV^, and under 5% elsewhere. To arrive to such a conclusion we considered the 
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dependencies on the angles 9k, and (p separately, i.e. in the way in which the three transverse 
asymmetries could be extracted if the B — ?• Kq£~^£~ events were absent. At large > 14 GeV^, 
instead, the effect of B ^ K^£+£- and Bs /o(980)£+£- on B ^ K*£+£- and Bs 4>£-^£- 
respectively, is completely negligible. 

This uncertainty, together with the one related to the charm loop [22], and the controllable 
uncertainties on the ratios of the B — )• K* form factors (R), suggests that the overall error on 
the transverse asymmetries (24) is under about 30%, and therefore at that level of accuracy 
their measurement remains a good tool for detecting the NP signal. This uncertainty can be 
circumvented only if the coefficient functions, appearing in the full distribution (18), are extracted 
from the simultaneous fit to the angular dependence of the decay events in three angles Oi, Ok 
and (j), which is a problem that could hopefully be solved experimentally once the number of 
events becomes sufficiently large for such a fit to be reliable. 
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A Explicit expressions 



In this appendix we collect the formulas that relate the coefficient functions xj^g (g^) appearing 
in the distributions (20,21,22) with the B — )• K* amplitudes ■M^^\Qt{<f) given in eq. (11), and 
with the B K^i+r amplitudes M^f \q^) given in eq. (13). We have: 



Xf(g2 



'o,t 

\\Mi\^ + \M[A^ + (L ^ Kf 



/32 . 



Am 



,2 . 



|MW|2 + 2Re(K'''K 



\Mi\^ + \M\\\^ + {L^ R) 



:^3(9') 

^4(g') 
2:5(9') 



2 

^/2 



|Xip-|A^{'|2 + (L^ii) 
Re (m^M\:*) +{L^R) 



V2pt [Re [Ml^Mi* - (L ^ fl))] , 



Re ( MI{M^* 



: \/2/3^ [im (m^M^* -{L^R) 
■ V2Pe Im (A^^A^f - {L ^ R) 

: ^ [Im (A^o^A^i*) + ^ ^)] 
: /3| [im (m^M^*) +{L^R) 



(31) 



Note that the above functions x\^'^\q^) are related to l\^''^\q^) of ref. [12] as z\^'''\q^) z[^'''\q^] 



(3/8) l\^'^\q^). Concerning the interference terms in eqs. (20,22) they read: 



/ A 

t ' 



•2 
J2 



(32) 
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